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Abstract 

The measurements of i?b = cr(e"'"e~ hh)/a{e'^e~ — > qq) and of the b quark 
forward-backward charge asymmetry, A^, at centre-of-mass energies above the Z 
pole are described. The measurement of -Rb is performed at between 130 and 
189 GeV using a b-tagging method that exploits the relatively large decay length 
of b-hadrons. The measurement of is performed using the large statistics event 
sample collected at a/s = 189 GeV with a lepton-tag analysis based on the selection 
of prompt muons and electrons. The results at -^i = 189 GeV are: 

Rb = 0.163 ± 0.013 (stal) ±0.005 (sysl), 
A^ = 0.61 ± 0.18 {stat.)± 0.09 (syst.). 
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Introduction 



The ratio R\, = a{e~^e~ bb)/cr(e+e^ qq) and the bottom quark forward-backward asym- 
metry are important parameters in precision studies of the Standard Model ^ and are 
sensitive probes to new physics. Their values have been measured very precisely at the Z 
pole PHOj- The measurements of these two quantities at centre-of-mass energies above the 
Z peak provide further tests of the Standard Model and put additional constraints on new 
physics, such as contact interactions [7]. 

In this paper the measurements of i?b at centre-of-mass energies between 130 and 189 GeV 
and of at y/s = 189 GeV obtained with the L3 detector [8^ at LEP are described. The 
main features that distinguish the production of bb pairs from lighter quark production are the 
long lifetime and hard fragmentation of b-flavoured hadrons and the large lepton momentum 
in semileptonic decays. The measurement of -Rb is based on a b-tagging method using lifetime 
information which uses data taken with the L3 Silicon Microvertex Detector (SMD) The 
measurement of A^ exploits the characteristic semileptonic decays of b-hadrons and relies on 
the good lepton identification and lepton energy resolution of the L3 detector. It requires 
leptons with high momentum along and transverse to the direction of the associated jet, caused 
by the hard fragmentation and high mass of the decaying hadron. Similar measurements have 
been published by other LEP collaborations |1II| . 

Event selection 

The data analysed in this paper are those collected by L3 from 1995 to 1998 at centre-of-mass 
energies between 130 GeV and 189 GeV. As the integrated luminosity at 130 GeV and 136 GeV is 
small, the corresponding data are combined and a luminosity weighted average centre-of-mass 
energy of 133.2 GeV is used. The centre-of-mass energies and the corresponding integrated 
luminosities used in the analysis are summarised in Table ^ 

The measurements of and A^ are performed using a sample of e"'"e~ Z*/Y — > qq(7) 
events selected with criteria similar to those used for the measurement of the e"'"e~ — > hadrons 
cross section ^1]. The signal events are required to have an effective centre-of-mass energy 
\/s' > 0.85a/s. This requirement rejects a large fraction of events with initial state radiation. 
Mis-reconstruction of the effective centre-of-mass energy induces a migration of radiative events 
to the kinematic region allowed by the cut on v^- This is taken into account in the analysis 
as an additional background, called radiative background. No correction is applied to the 
final results to take into account the interference between initial- and final-state radiative 
corrections. 

To further reject radiative events, the visible energy, ii^vis, is required to be greater than 
0.6^/s and the longitudinal energy imbalance must be less than 0.25i?vis- The events are 
selected in a fiducial region defined by | cos6 \ < 0.85, where 6 is the polar angle of the thrust 
axis. This cut ensures good track quality and efficient lepton identification. For centre-of-mass 
energies larger than the W-pair production threshold, additional cuts are needed to reject the 
W"'"W~ background. To reject hadronic W decays, the variable is required to be less than 
0.011, where 1/34 is the Durham jet finding algorithm parameter ^2] for which the transition 
from three to four jets occurs. As the W~'"W~ background increases for higher centre-of-mass 
energies, the maximum allowed value of 7/34 is lowered to 0.008 at y/s = 183 GeV and to 0.006 
at 189 GeV. Furthermore, to reject semileptonic W decays, the transverse energy imbalance 
must be smaller than 0.4i?vis- For the A^ analysis the presence of an identified electron or 
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muon with an energy smaller than 40 GeV is required. 

Efficiency and background studies are performed for each centre-of-mass energy using the 
following event generators: PYTHIA T3j (e"'"e~ — > hadrons(7), Zee), KORALZ 14^ (r+r^(7)), 
PHOJET (hadronic two-photon collisions), KORALW [10] (e+e" ^ W+W-(7)) and EX- 
CALIBUR [T7] (e+e~ — > qqq'q'), which includes the ZZ production diagrams. 

Measurement of 

The measurement of i?b is based on a b-tagging algorithm that exploits the relatively large 
decay length of b-hadrons JH]- The confidence level, C^, that a set of tracks originated 
from the primary vertex, is constructed using the decay length significance L/ai of each track. 
First the crossing point of each track with the closest jet is determined in both the r0 and 
sz projections. Then the signed distances {dnj,, dgz) between these crossing points and the 
reconstructed primary event vertex are projected onto the jet axis to determine the decay 
length, L. If the probability that both r0 and sz measurements are compatible exceeds 5%, 
then the two are combined. Otherwise, only the r0 projection is used. 

Small biases in drtp and dgz are removed by recalibrating their mean values as functions 
of the azimuthal and polar angles of the tracks. The errors on drtp and dsz are parametrised 
to take into account the contributions coming from the track fit, the primary vertex position 
measurement and the effect of multiple scattering. The estimated errors on dr^ and dgz are 
rescaled by comparing their values with those observed in a sample of tracks with negative decay 
length. The scale factors are found to be within 10% of unity. The tracking resolution depends 
critically on the pattern of SMD hits associated to the tracks. Therefore the tracks are grouped 
into different classes according to the vertex detector hit pattern and for each class different 
resolution parameters are considered. The study of the tracking resolution described above is 
performed using samples of hadronic Z decays corresponding to an integrated luminosity of 
about 2.5 pb~^ each, collected each year during calibration runs at the Z peak. 

The confidence level, Cn, is calculated by taking into account the fraction of tracks with 
positive decay length. 



where A^+ is the number of tracks with positive decay length. The probability that a track 
originated from the primary vertex, P(L/(Tl), is obtained by fitting the distribution of the 
decay length significance for tracks with negative decay length. The shape of this distribution 
is due to detector resolution effects and is modelled by a resolution function parametrised as 
the sum of two gaussians and an exponential tail. The distributions of the discriminant variable 
D = — logiQ(CAr) obtained at the different centre-of-mass energies are shown in Figure [U 

The measurement of -Rb is performed using an event-tag method by counting the events 
containing b quarks. The b-events are selected by applying a cut on the discriminant where the 
cut position is chosen so that the expected statistical error on the cross section cr(e"'"e~ — > bb) 
is minimised. Given the number A^"*^^ of tagged events and N^^^ of tagged background events, 
i?b is derived from the formula: 




(1) 



1 [ATf ^ - N^^^ 




(2) 
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where N and N ^ are respectively the total number of selected qq events in data and the 
total number of background events from Monte Carlo; e^, Ec and ^uds are the tagging efficiencies 
for b, c and light quarks, respectively, and are estimated from Monte Carlo. 

The value of the cc production cross section, R^, is taken from the Standard Model prediction 
calculated by ZFITTER ^19j. The dependence of -Rb on can be parametrised as 

Ai?b = a{R^)AR^ (3) 

where ARc = Rc — Rf^. The values of Rf^ used in the measurement and the corresponding 
coefficients a(-Rc) are shown in Table |21 

The background includes the contribution from non-qq events and the contamination of 
radiative hadronic events. The radiative background decreases from 16% at 133 GeV to 4% at 
189 GeV. The background due to processes other than e"'"e~ qq(7) is dominated by hadronic 
W decays for > 161 GeV and varies from 1% at = 161 GeV to 11% at 189 GeV. The 
background from e"'"e~ — qqq'q' events is at the level of 1% at 189 GeV and negligible at lower 
centre-of-mass energies. The contribution from other processes is negligible. 

The number of selected and tagged events and the background contamination both in the 
total and in the tagged sample for each centre-of-mass energy are shown in Table |H1 The 
tagging efficiencies are also quoted along with their statistical errors. The b tagging efficiency 
decreases slightly with increasing centre-of-mass energy due to the larger non-qq background 
contribution. The measured values of i?b are shown in Table El together with their statistical 
errors. 

The stability of the measurements has been checked by varying the cut on the b-tagging 
discriminant. The results of this check for the two high statistics samples at ^/s = 183 GeV 
and 189 GeV are shown in Figure El No statistically significant deviation from the i?b value 
obtained with the nominal cut is observed. The event-tag method used in this analysis has 
been cross-checked by measuring i?b using the calibration data collected each year at the Z 
peak. The results are found to be in good agreement with our published value jlj. 

Systematic uncertainties on 

The largest contribution to the systematic uncertainty on i?b comes from the description of 
the tracking resolution function in the Monte Carlo, since it affects the estimation of the 
tagging efficiencies involved in the i?b measurement. In order to study this source of systematic 
uncertainty, several Monte Carlo samples were reconstructed, where the parameters used to 
describe the resolution function have been varied according to their estimated uncertainty. The 
differences between the values of ^b, £c and ^uds obtained from these Monte Carlo samples and 
the nominal values shown in Table have been propagated to obtain an estimation of the 
systematic uncertainty due to tracking resolution effects. 

Another source of systematic uncertainty comes from the description in the Monte Carlo 
of the tracking efficiency and of the relative fraction of tracks with a given SMD hit pattern. 
To study this effect, the tracking efficiency has been varied in the Monte Carlo by ±2% and a 
migration of tracks at the level of 1% between the different classes has been introduced. The 
systematic uncertainty coming from this source has been obtained as above by propagating the 
variations on the flavour tagging efficiencies. The systematic uncertainties discussed above are 
added in quadrature to give the total systematic uncertainty due to tracking effects. 

The modelling of b and c quark fragmentation and decays introduces a systematic uncer- 
tainty on i?b since the tagging efficiency for b and c quarks is obtained from Monte Carlo. 
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Modelling uncertainties of the b-hadron properties are estimated by varying the mean value 
of the b-hadron energy fraction (xeO^)), the charged decay multiplicity and the average life- 
time of b-hadrons. For the estimation of Ec an accurate knowledge of production and decay 
properties of the c-hadrons is important, since the different species, D°, D"*", Dg and A^., have 
lifetimes varying in the range of 0.2 to 1.1 ps ji20j. The variation of the parameters is performed 
following the recommendations of Reference [21] . A breakdown of the systematic uncertainties 
due to the modelling of b and c-hadron physics is given in Table IH 

The error due to the finite Monte Carlo statistics has been estimated by propagating the 
statistical errors on eb^^c and e^ds- The systematic uncertainty due to the selection of the 
e"'"e~ — > qq(7) event sample has been estimated by varying the selection cuts and taking the 
corresponding variation of i?b as a systematic uncertainty. The error coming from the non-qq 
and radiative background estimation was evaluated by varying the background fractions within 
their estimated uncertainties. 

The various contributions to the systematic uncertainty on i?b are summarised in Table El 

Measurement of 

The forward-backward charge asymmetry of b quarks is measured by tagging e"'"e~ — >■ bb events 
with electrons and muons with high momentum and high transverse momentum with respect 
to the nearest jet. The leptons are also used to identify the charge of the b quarks. The quark 
direction is estimated using the event thrust axis. To enhance the b purity, the b-tagging 
method previously described is used. 

Lepton identification and event selection 

Muons are identified in the muon chamber system. The muon tracks must include track seg- 
ments in at least two of the three layers of muon chambers and must point to the interaction 
region. 

Electrons are selected in the barrel region | cos^^ | < 0.7. They are identified by an energy 
deposit in the electromagnetic calorimeter consistent with an electromagnetic shower matching 
with a track in the central tracker. To reject misidentified hadrons, the energy in the hadron 
calorimeter behind the electromagnetic shower is required to be less than 6 GeV. 

The momentum of muon candidates is required to be greater than 6 GeV, while that of 
the electron candidates must be greater than 3 GeV. The lepton transverse momentum, pt, is 
defined with respect to the nearest jet, where the measured momentum of the lepton is not 
included in the jet reconstruction. If there is no jet with an energy greater than 6 GeV in the 
same hemisphere as the lepton, the pt is calculated relative to the thrust axis of the event. The 
transverse momentum of the leptons must be greater than 0.7 GeV to enhance the bb purity. In 
events containing more than one lepton candidate, only the lepton with the highest pt is used 
in the asymmetry analysis. In order to further reduce the background coming from radiative 
qq events and W pair production, the angle a between the lepton and the event thrust axis 
must satisfy the requirement | cos a | > 0.95. 

A cut in the {pt, D) plane is applied to enhance the b purity, where D is the discriminant 
variable used for the measurement of i?b- The cut is defined by: 

D > -1.8pt + 2.4 (muons), 

D > -0.7pt + l-6 (electrons). (4) 
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These cuts were chosen as a resuh of an optimisation procedure. 

The number of prompt lepton candidates is 136, comprising 71 muon candidates and 65 
electron candidates in good agreement with the Monte Carlo expectation of 139.4 events with 
75.3 muons and 64.1 electrons. 

Monte Carlo events are classified into eleven categories as listed in Table El Leptons in b 
events not arising from the semi-leptonic decays of categories 1 to 5 are classed as "b fake 
lepton" . This set of events consists of leptons from vr and K decays, Dalitz decays and photon 
conversions and of misidentified hadrons. The Monte Carlo estimate of the sample composition 
is shown in Table IHl together with the asymmetry contribution of each event class. 

The purities are 33.2% for the muon sample and 16.8% for the electron sample. The effi- 
ciencies corresponding to these purities are 21.2% and 7.6% respectively. The momentum and 
transverse momentum spectra for the selected muon and electron candidates together with the 
sample composition are shown in Figure El 

Asymmetry determination 

The b quark scattering angle is obtained from the observable cos^b = —Q cos6't, where the 
thrust direction 6^ is oriented into the hemisphere containing the lepton and q is the lepton 
charge. The distribution of cos^b for the selected muon and electron samples is shown in 
Figure jH 

Monte Carlo studies of the "b fake lepton" class indicate a residual correlation between 
the quark charge and the observed charge, for both muons and electrons. Variations of the 
generator level asymmetry reveal a linear correlation with the observed asymmetry for this 
class of events with y4bckg=(0.16 ± 0.09) xAb for muons and (0.25 ± 0.08) xAb for electrons, 
where the errors on the factors are statistical. Variations by the statistical errors on the factors 
are used to estimate systematic uncertainties arising from the Abckg assignment. A comparison 
with results obtained using a constant Abckg indicate differences of 0.013 and 0.018 in the final 
measured asymmetries for muons and electrons, respectively. 

For each lepton i in the data, the probability /fc(«) to belong to each of the first six categories 
listed in Table (HI is determined from the number and type of Monte Carlo leptons found in the 
appropriate cos6'b bin. The asymmetry is obtained by applying a maximum-likelihood fit 
assummg no mixing. The likelihood function has the form: 

. _ n^r eLi + cos^ o^) + Ak cose,] 

nft^ (ELi + cos2 e,) + A, cos^,])^0) 

where A^ is the asymmetry contribution for the signal category k and W{j) represents the 
Monte Carlo to data normalisation relevant to the Monte Carlo lepton j, obtained from the 
ratio of the total number of data and Monte Carlo leptons found in the corresponding cos^b 
bin. The product over non-b Monte Carlo leptons in the denominator takes into account non-b 
background contributions in the data sample. 

The fitted asymmetries, corrected for the charge confusion (2.0 ± 0.2% for electrons, and 
0.1 ± 0.1% for muons) and B°B° mixing (xb = 0.1192 ± 0.0068 |11E2|), are 

= 0.70 ±0.22 (stal), 
A)^^ = 0.39 ± 0.34 (stal). (6) 
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Systematic uncertainties on A 



The dependence of the measured asymmetry on other electroweak parameters -Rb, Rc or 
used in the Monte Carlo modelling can be parametrised by : 

AA^ = a{X)AX (7) 

where AX = X - X^^ for X = i?b, or A^. The values of R^^, i?^M ^c,sm ^^^^ ^-^^ 
Monte Carlo are given in Table [7| together with the values of the coefficients a{X). 

Relevant fragmentation parameters, semileptonic decay models, and branching ratios are set 
to their measured values and their uncertainties are used to estimate the systematic uncertainty 
on the asymmetry. The parameter values and variations follow the recommendations developed 
in Reference |2I] and are listed in Table |H1 along with the corresponding systematic uncertainties 
on <. 

Effects due to uncertainties in the lepton momentum resolution are estimated by smearing 
the muon momentum by 1% and the electron momentum by 5%. The uncertainty on the charge 
confusion correction is included as a contribution to the systematic uncertainty. 

A possible bias on the measured asymmetry due to the non-b background contribution is 
estimated by varying each background fraction by ±5% and the light quark asymmetry by 
±0.02. The corresponding variations in A^ are taken as systematic uncertainties. 

To study the Monte Carlo description of the pt spectrum of the non-b background, Monte 
Carlo events have been compared with high statistics background-enhanced data samples. 
Using these data samples to constrain the shape of the pt spectrum, Monte Carlo reweighting 
factors for non-b background are obtained in each pt bin and used to modify the weights W{j) 
in Equation |S| Half of the change in the measured asymmetry caused by this reweighting 
procedure is assigned as a systematic uncertainty from the background pt description. 

The total systematic uncertainty on A^ is given in Table |H1 

Results 

The measured values of i?b are summarised in Table El with their statistical and systematic 
errors. The measurements are compared to the Standard Model prediction in Figure At 
= 189 GeV the measured value of i?b is: 

Rb = 0.163 ± 0.013 (stat.) ± 0.005 (syst.) 

to be compared with the Standard Model value [T^ of 0.166. 

The two measured values for A^'^ and A^^ (EquationlHI) are consistent and they are combined 
to give at y/s = 189 GeV: 

A^ = 0.61 ± 0.18 (stat.) ± 0.09 (syst.) 
to be compared with the Standard Model value [ini of 0.58. 
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(GeV) 


133.2 


161.3 


172.1 


182.7 


188.6 


L (pb-i) 


12.01 


10.95 


10.25 


55.65 


176.3 



Table 1: Summary of centre-of-mass energies and integrated luminosities considered 
in this analysis. 



x/i (GeV) 


133 


161 


172 


183 


189 


rip 


0.221 


0.242 


0.246 


0.250 


0.252 




-0.113 


-0.118 


-0.106 


-0.099 


-0.100 



Table 2: Values of i?^'^ used for the measurement of i?b and coefficients a(i?c) giving 
the dependence of the measured i?b on (Equation IH}. 



(GeV) 


jyobs 


jyobs 








tagging 
efficiencies 














^b 


=0.676±0.007 




133 


850 


134 


138.5 ±2.0 


23.8 ±0.4 


^uds 


=0.097±0.004 
=0.023±0.001 


0.177±0.023 












£b 


=0.645±0.006 




161 


319 


45 


28.2 ±0.4 


3.9 ±0.2 


^uds 


=0.103±0.003 
=0.031±0.001 


0.152±0.035 












^b 


=0.619±0.006 




172 


266 


43 


38.0 ±0.4 


3.7±0.1 




=0.095±0.003 
=0.032±0.001 


0.212±0.045 












£^b 


=0.621±0.004 




183 


1172 


147 


185.8 ±0.9 


16.6 ±0.3 


£c 
^uds 


=0.091±0.002 
=0.033±0.001 


0.145±0.020 












£b 


=0.600±0.002 




189 


3462 


486 


524.3 ±2.3 


50.6 ±0.6 


£c 
^uds 


=0.099±0.001 
=0.044±0.001 


0.163±0.013 



Table 3: Number of selected (A^°''^) and tagged (A^"''^) events, background con- 
tamination both in the total {N'°'^^^ and in the tagged sample (A^]?''^^) and tagging 
efficiencies for each centre-of-mass energy. The errors on the tagging efficiencies are 
statistical. In the last column the measured values of i?b are listed. The errors are 
statistical only. 



13 



Source 




133 G(¥ 


161 Gc¥ 


172 GeV 


183 GeV 


189 GcV 


b fragmentation 


0.0015 


0.0014 


0.0020 


0.0014 


0.0016 


b lifetime 


0.0003 


0.0003 


0.0004 


0.0003 


0.0003 


b decay multiplicity 


0.0010 


0.0009 


0.0013 


0.0009 


0.0011 


c modelling 


0.0005 


0.0006 


0.0006 


0.0006 


0.0006 


Total b, c physics modelling 


0.0019 


0.0018 


0.0025 


0.0018 


0.0021 



Table 4: Summary of systematic uncertainties on i?b due to b and c-hadron physics 
modelling. 



Source 




133 GeV 


161 GeV 


172 GeV 


183 GeV 


189 GeV 


Tracking effects 


0.004 


0.004 


0.005 


0.003 


0.004 


b, c physics modelling 


0.002 


0.002 


0.003 


0.002 


0.002 


MC statistics 


0.003 


0.002 


0.003 


0.001 


0.001 


Event selection 


0.002 


0.002 


0.002 


0.002 


0.002 


Background 


0.001 


0.001 


0.001 


0.001 


0.001 


Total systematic uncertainty 


0.005 


0.005 


0.007 


0.004 


0.005 



Table 5: Summary of systematic uncertainties on R^,. 
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Category 


Muons (%) 


Electrons (%) 




b^l 


33.20 ± 0.40 


16.80 ± 0.30 






8.20 ± 0.20 


1.45 ± 0.10 


-A 


1) ^ c ^ 1 


1.80 ± 0.10 


0.15 ± 0.05 


A, 


b ^ J/^/- ^ 1 


0.40 ± 0.05 







b ^ r ^ 1 


1.10 ± 0.09 


0.70 ± 0.06 




b — s> fake 1 


5.70 ± 0.20 


12.90 ± 0.30 


^bckg 


c^l 


14.70 ± 0.35 


2.10 ± 0.10 




c fake 1 


5.50 ± 0.40 


12.70 ± 0.45 




u, d, s — > 1, fake 1 


21.10 ± 0.25 


37.20 ± 0.30 




rad. background 


3.80 ± 0.20 


4.40 ± 0.10 




WW background 


4.50 ± 0.20 


11.30 ± 0.45 





Table 6: Fractions of events in the various categories described in the text, esti- 
mated from Monte Carlo. The corresponding asymmetry contributions used in the 
likelihood fit are shown in the last column. 
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Standard Model value 


a{X) 


R^"^ = 0.166 
i?^M = 0.252 

A"^^"^ = 0.62 


-0.81 
+1.44 
+0.47 



Table 7: Values of the Standard Model parameters R^,^, Rf^ and used in 

the analysis and coefficients a(X) giving the dependence of A^ on them (Equa- 
tion I^). 



Contribution 


Value 


Variation 




XE{h) 


0.702 


± 0.008 


+ 0.004 


xe{c) 


0.484 


± 0.008 


± 0.006 


Br(b ^ 1) 


0.105 


± 0.005 


+ 0.007 


Br(b c ^ 1) 


0.080 


± 0.005 


± 0.004 


Br(b ^ c ^ 1) 


0.013 


± 0.005 


+ 0.001 


Di\p — ^ r — > Ij 


O.OOo 


1 n nn 1 
± O.OOi 


-T- n nno 
+ 0.002 


Br(c ^ 1) 


0.098 


± 0.005 


± 0.013 


Br(b ^ J/^A ^ 1) 


0.001 


± 0.001 


± 0.002 


Fragmentation and Branching 


Ratios 




0.017 


b 1 model 






+ 0.005 


c 1 model 






± 0.025 


b ^ D model 






± 0.003 


Decay Models 






0.026 


background fractions 




±5% 


± 0.024 


u, d, s background asymmetry 


0.00 


±0.02 


+ 0.020 


^bckg asymmetry 




±0.08 


+ 0.004 


Background effects 






0.031 


Pt background description 






± 0.073 


lepton momentum smearing 






± 0.005 


charge confusion correction 






+ 0.002 


mixing dependence 


0.1192 


0.0068 


± 0.010 


Detector effects 






0.074 


Total systematic uncertainty 






0.087 



Table 8: Systematic uncertainties on Al^. The values and variations for the charge 
confusion correction and momentum smearing are given in the text. 
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(GeV) 


Rh 


SM Expectation 


133.2 


0.177 ±0.023 ±0.005 


0.185 


161.3 


0.152 ±0.035 ±0.005 


0.172 


172.1 


0.212 ±0.045 ±0.007 


0.169 


182.7 


0.145 ±0.020 ±0.004 


0.167 


188.6 


0.163 ±0.013 ±0.005 


0.166 



Table 9: Summary of results on R^. The first error is statistical, and the second is 
systematic. The values predicted by the Standard Model [11] are shown in the last 
column. 
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Figure 1: Distributions of the b-tagging discriminant for the selected events at 
the different centre-of-mass energies. The points are the data and the sohd hnes 
represent the total Monte Carlo expectation. The open histograms show the sig- 
nal contribution, the vertically hatched histograms represent the c and light quark 
background from qq events and the cross-hatched histograms represent the sum of 
non-qq and radiative background, as indicated in the legend of the upper left figure. 
In the bottom right plot, the distribution for the total selected sample is shown. 
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Figure 2: Variation of Ri, as a function of the cut on the b-tagging discriminant for 
the high statistics data samples at = 183 GeV and 189 GeV. The sohd circles 
correspond to the measured values of R^. 
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Figure 3: Momentum and transverse momentum distributions for selected leptons. 
The histograms show the estimated composition of the sample according to Table IHl 
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Figure 4: Distribution of the scattering angle cos^b = — qcos^T for the muon and 
electron samples, respectively. The histograms represent the sample composition 
according to Table El 
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Figure 5: Measured values of i?b as a function of centre-of-mass energy. The curve 
is the Standard Model prediction [T^ . 
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